The potential impact of groove modes on Type II planetary 

migration 
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ABSTRACT 

In this letter, we briefly describe the evolution of a variety of self-gravitating protoplanetary 
disk models that contain annular grooves (e.g. gaps) in their surface density. These grooves 
are inspired by the density gaps that are presumed to open in response to the formation of a 
giant planet. Our work provides an extension of the previously studied groove modes that are 
known in the context of stellar disks. The emergence of spiral gravitational instabilities (GI) is 
predicted via a generalized eigenvalue code that performs a linear analysis, and confirmed with 
hydrodynamical simulations. We find the presence of a groove drives a fast-growing two-armed 
mode in moderately massive disks, and extends the importance of self-gravitating instabilities 
down to lower disk masses than for which they would otherwise occur. We discuss the potential 
importance of this instability in the context of planet formation, e.g. the modification of the 
torques driving Type II migration. 

Subject headings: stars: planetary systems - hydrodynamics 



1. Introduction 

The theory of giant planet formation has 
evolved very rapidly over the past decade, and 
the now near-p aradigmatic core accretion theory 
(e.g. iPollack et al.lll996l ) has been greatly refined 
and endowed with a var i ety of physical inputs (see, 



relati on (pantos et al.l l2003t iFischer fc Valenti 
2005). Another criticism is that GI requires the 



e.g. lAlibert et alj|2004 iHubickvi et al.ll2005fh It 



is now generally accepted that Jupiter-mass plan- 
ets grow by initially accumulating 5-10 M^ cores, 
which then accrete large quantities of gas from the 
surrounding nebula. During the rapid gas accre- 
tion phase, the growing Jovian planet opens a gap 
in the disk and its orbit subseque ntly evolves via 



the process of Type II migration (|Lin et al 
Papaloizou et alj|2007h . 



19961 : 



A competing theory (e.g. lBossll2000l ) holds that 
giant planets such as Jupiter fragmented directly 
from the nebula as a result of gravitational insta- 
bilities. This mechanism may account for some of 
the observed extrasolar planets, but it has a num- 
ber of difficulties if invoked to account for the full 
population. In particular, it fails to explain the 
observationally well established planet-metallicity 



Toomre Q parameter, Q — Kc s /irGa, to be of 
order unity at some radius in the disk. This, 
however, generally requires rather massive disks, 
Mr)/M s tar j^ 0.1, which are not often observed. 

Global GIs depend not just on disk mass and 
sound speed, but al so on the densit y profile of the 
disk. As shown by iToomrd ( 1981 ). a sharp den- 
sity gradient (i.e. "edges") can drive a wave cy- 
cle that leads to th e subsequent destructio n of the 
edge. Analogously, ISellwood fe~Lml (jl989h [SL89] 
described a family of GIs in stellar disks dubbed 
groove modes which comprise fast-growing distur- 
bances driven by narrowly defined structures in 
particle angle-action space. These features, in 
physical space, correspond to narrow density de- 
pressions; groove modes are thus akin to the edge 
mode in that they are driven by steep density 
gradients at corotation. SL89 further envisioned 
a feedback cycle which kicks in when a slowly- 
growing intrinsic mode (i.e. an inner edge mode) 
carves a groove by scattering particles in a nar- 
row range at its Lindblad resonance and creating 



a fast-growing groove mode; the new groove gives 
rise to further groove modes once it travels to its 
own Lindblad resonances. 

This feedback cycle relies on the wave-particle 
interaction at the resonances, and thus cannot 
be immediately realized in a gaseous disk. The 
groove mode may have relevance, however, in the 
protoplanetary con t ext: t h eoretical considerations 
(ILin fc Papaloizoul 1979t Goldreich fe Tremaind 



1980c ILin fe Papaloizoul 19931). numerical simu 
lations (e.g. iTakeuchi et al.l 11996c iBrvden et al. 

1999) and tentative observational indications 



(jSetiawan et al.ll2008f ) point to the fact that mas- 
sive planets can resonantly drive trailing waves 
that transport angular momentum and open a 
"gap" in the disk. The detailed process of open- 
ing and maintaining the gap depends on the bal- 
ance between the angular momentum flux result- 
ing from spiral waves driven by the planet and 
that due to the viscosity of the disk. When the 
planet is sufficiently massive and the viscosity 
is low enough, a stable groove can be naturally 
carved and mantained on the disk. 

In this letter, we open an investigation into 
the possibility that the planet-induced density 
gap can drive a fast-growing, self-gravity induced 
groove mode in a disk, sidestepping the prob- 
lems of carving a groove in the first place and 
providing a mechanism for maintaining the am- 
plifier. We support our hyp othesis with both a 
linear mode analysis (see, e.g. I Adams et al.l 1 19891 ) 
and a full hydrodynamical simulation (following 
Laughlin fc Rozvczkalll996l hereafter LR96) of a 
disk with an imposed groove in surface density, 
as inspired by prev ious numerical inv estigations 
for stellar disks f e.g. ISellwood fe Kahn|[l99ll ). We 
span a range of mass ratios to assess the relevance 
of this GI to the realistic disk masses observed. 

2. Procedure 

We employ a two-dimensional hydrodynamical 
grid code for following the evolution of a thin, 
self-gravitating disk. The continuity and Eu- 
ler equations in polar coordinates are solved us- 
ing a second-order van Leer type scheme, cou- 
pled with time stepping that is first-order accu- 
rate. The basic differen ce equations are given in 
Stone fc Normanl (jl992h . The self-gravity of the 
disk is obtained by applying the Fourier convo- 



lution theorem t o the potential dictated b y the 
Poisson equation ( Binnev fc Tremainelll987l ). The 
de tails of the hyd r odyn amical code are described 
in lLaughlin et al.l (J1997P and related papers. 

We adopt the following parametrization for the 
surface density of the disk: 



a(r) = a e^ r - R ^n 



AA 2 



Rp) 



which represents a Gaussian profile multiplied by 
a Lorentzian gap of depth A, semi-width A and 
central position Ftp; for A = 0, this profile is 
the "reference disk" considered by LR96. The 
disk model used throughout this letter has no 
pretense of actually representing a protoplanetary 
disk faithfully, but has the advantage of possess- 
ing a single intrinsic m = 2 mode for A = 
that is clearly identifiable both in semianalytic cal- 
culations and in the nonlinear simulations. The 
choice of a Lorentzian profile to represent the gap 
i s arb itrary as well, and follows ISellwood fc Kahn 
(|l99ll ). We take w = 0.03, R = 0.25 and an inner 
edge of 0.05. A polytropic equation of state is as- 
sumed, P = Ka 1 with 7 = 2 (again, this is chosen 
largely for illustrative purposes). The character- 
istic width A is meant to represent a typical gap 
width; the chosen value of 0.07 can be derived from 
the approx imate scaling derived bv lVarniere et al. 
(e.g. l2004h with n ~ 5 x 10 5 and q = 2 x 10" 3 . 



The set of units used in the code takes the outer 
grid radius Rd and the gravitational constant G 
equal to unity and M* = 0.5. 

To quantify the strength, pattern speed and 
growth rate of each spiral mode, we compute the 
Fourier decomposition of the surface density, de- 
fined as 



1 
27 



(7(r,$)e-"" ? d$ = Cm x ffl o, (2) 



for a mode number m. A global measure of the 
growth of a particular mode is given by integrat- 
ing the m-th Fourier amplitude a m over the radial 
range and normalizing to the azimuthal average of 
the surface density: 



In 



pRd 

JRi 



i(r)dr 



L D a {r)dr 



elf 



log C„ 



(3) 



The phase of a disturbance is recovered as 



$ m (r) = tan l [Im(— a m )/Re(a m )]; the local pat- 
tern speed is then given by Sip = (l/m)$ m . 

The growth rate and pattern speed that emerge 
from the hydrodynamical simulation are checked 
against a linear numerical analysis code w e devel- 
oped, as described e. g. in lLaughlin et al.1 (|1997l ). 
which solves a matrix equation akin to a gener- 
alized eigenvalue problem; the solution is valid in 
the linear regime and yields a complex eigenvalue, 
which indicates the pattern speed Sip and growth 
rate 7, and a complex eigenvector, which describes 
the radial variation and phase of the mode. A 
graphical depiction of the procedure is shown in 
Figure ([I]). Comparison with the full nonlinear 
simulations enables us to check the consistency 
and accuracy of the two independent approaches. 

3. Computer simulations 

Table ([1]) lists the disk models considered in 
this letter. We have first set up a "base" (A = 0) 
disk in rotational equilibrium with surface density 
given by Equation [I] with qr> = Mn/M s t a r = 1 
and Qmin =1.21 (Model 1); this sets the normal- 
ization constant <7o and the polytropic constant 
K. The resulting Toomre Q profiles are shown in 
Figure ^. The equilibrium was disrupted with 
a random density perturbation of order 10 _3 cr(r). 
The grid covers the polar coordinates (r, <p) with 
256 logarithmically spaced zones in radius and 256 
equally spaced azimuthal zones. Each model is 
evolved for at least 100 time units, although the 
reflective boundaries muddle the nonlinear evolu- 
tion once the wave reaches the outer part of the 
disk. 

Figure ([3]) shows the evolution of the sur- 
face density normalized to the azimuthal aver- 
age £(r) = log(a(r)/a(r)). In accordance with 
the previous investigations and with our linear 
code, the base disk is unstable to a single m = 2 
"grand-design" spiral mode. The GI grows in the 
linear regime for the first few dynamical times, 
and within about 10 dynamical times it visibly 
perturbs the outer edges of the disk. After a 
few more dynamical times, the spiral pattern has 
reflected off the boundaries and has propagated 
back into the densest regions of the disk. The nor- 
malized amplitude C2 (Figure Q| shows that the 
dominant two-armed mode grows at an exponen- 
tial rate (linear regime) until stabilizing around 



a constant amplitude (mode saturation). For our 
purposes, the linear growth rate is the primary 
quantity of interest. 

A second set of disk parameters (Model 2) in- 
cluding a density groove was set up, with A — 0.90 
(a 90% dip), A = 0.07 and R P = 0.4. Maintaining 
the same normalizations for density and pressure 
as above, the disk has the same density profile suf- 
ficiently far from the groove, but a smaller disk- 
to-star mass ratio, qo = 0.63. By visual inspec- 
tion of Figure ([3]) and the slope of the normalized 
Fourier amplitude in Figure ([4]), it is clear that 
the presence of the groove drives a far more vio- 
lently growing instability. The dashed linear slope 
is derived from the mode analysis for the same 
disk parameters, and agrees well with the hydro- 
dynamical code; the small discrepancies between 
the two methods can be traced to the effective soft- 
ening given by the FFT-based solution scheme for 
the Poisson equation. We then varied the density 
normalization to yield qo = 0.32, 0.16, 0.13, 0.08, 
and 0.06 (Models 3-11). 

Measuring the pattern speed is more complex, 
since the spiral rotates only approximately in a 
rigid fashion. We thus calculate a "global" mea- 
sure Sip by weighing the local Slp(r) with the lo- 
cal density enhancement |c m |, so that the densest 
parts of the spiral contribute the most to the pat- 
tern speed. The corotation radius of the models 
with a groove is found to lie within the gap, con- 
firming its nature as a groove mode and in accor - 
dance to what found bv lSellwood fc Kahnl (J1991I ). 

Since the presence of the groove diminishes the 
qu parameter for a given density normalization, 
we also set up a series of base disks (Models 12- 
17) with a density normalization chosen so that its 
mass would equal that of Models 2-7. For a given 
qo , in the presence of the groove the m = 2 mode 
exponentiates about 3-4 times faster than the base 
disk. Models 14-17, despite having equal disk mass 
and similar Q m in than their groovy counterparts, 
do not show an exponential phase during the du- 
ration of the simulation. A growing groove mode 
is detectable down to qr> = 0.08 (Model 6), but 
does not show a resolvable exponential phase for 
q D = 0.06 (Model 7). The m = 2 groove instabil- 
ity in Model 6 is growing about as fast as Model 
13, which is four times as massive. 

The natural outcome for the unstable models 
is the transport of mass and angular momentum, 



resulting in the filling of the gap. The evolution of 
the azimuthally averaged density profile for the 
two lowest-mass unstable models (Model 5 and 
6) is shown in Figure $$. An approximate esti- 
mate for the effective ct-typ e viscosity coefficient 
(jShakura fc Svunvaevlll973l ) is derived by a proce- 
dure similar to that employed in LR96, although 
the effective viscosity given by GIs is not well 
characterized by a local prescription. We solved 
the diffusion-type equation for the time-dependent 
surface density and compared it with the evolution 
in the hydrodynamical simulation, once the spiral 
pattern has established itself. This simple estima- 
tion yields a ~ 0.16 for qp = 0.13 and a ~ 0.04 
for qo — 0.08, and agrees approximately with the 
timescale for closing the gap via viscous diffusion. 

4. Discussion and conclusion 

In this letter, we have found that disk self- 
gravity may play a significant and as-yet largely 
unstudied role in disks in which a planet has 
opened a gap in the surface density profile of the 
disk. Gravitational torques from massive proto- 
planets necessarily impart a surface density gap 
in a disk, and the formation of a gap provides 
the struc ture needed b oth for a feedback ampli- 
fier (e.g. lToomrelll981 ) as well as for the groove 
mode (SL89). In essence, a gap provides a pres- 
sure gradient which can locally reduce the effec- 
tiveness of the epicyclic frequency in stabilizing 
the disk against its own self gravity, thus allowing 
instability at low surface densities. In the absence 
of a perturbing planet, the nonlinear outcome of 
a groove instability would be to destroy the sharp 
density gradient that promoted the instability in 
the first place. In a planet-forming disk, however, 
the embedded Jovian planet will exert torques 
whose net effect is to maintain the gap. The com- 
petition between the Gl-induced torques and the 
planetary torques may thus lead to a significant 
modification of the criterion for gap opening (as 
evidenced by the large effective viscosity found in 
our experiments) , which in turn can have a signif- 
icant effect on the resulting migration and growth 
of the protoplanet. We have found that this effect 
is likely to be relevant even when the disk mass is 
lower than the qr, > 0.1 — 0.2 value at which signif- 
icant GIs are g enerally thought to occur (see, e.g. 
Shu et all [19901 JBossI 1l997l . LR96). By compari- 
son, the typical mass within 30 AU assumed for 



the Minimum Mass Solar Nebula yields qu ~ 0.01 
(this likely underestimates by a f a ctor o f 3) . Fur- 
ther afield, lAndrews fc Williams! (|2005h surveyed 
153 young stellar objects in the Taurus Aurigae 
star-forming complex, and found a median disk 
mass of qo = 0.5%. While these disks span a va- 
riety of ages, and show a variety of masses, it is 
not clear how massive the average protoplanetary 
disks is at the time when cores enter the rapid gas 
accretion phase, and whether the mass might be 
enough to trigger the instability studied here. It 
must also be kept in mind that while the idealized 
disks in the simulations shown here display no GI 
for qo < 0.08, the limited numerical resolution 
and high intrinsic numerical viscosity allows us to 
identify only modes with intrinsically rapid growth 
rates. Further work with realistic disk models is 
required to find the true minimum mass for the in- 
stability. We are eager to continue this analysis by 

(1) including the planetary potential in both our 
linear analysis and in our nonlinear simulations, 

(2) carrying out the simulations at higher numeri- 
cal resolution, and (3) adopting more realistic disk 
models. 
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Origins of Solar Systems Program through grant 
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Table 1 
The table lists growth rate (ju ni %i) and pattern speed (Ct P , Cl P ) as measured 

RESPECTIVELY BY THE LINEAR CODE AND THE FULL HYDRODYNAMICAL SIMULATION FOR THE VARIOUS 
DISK MASSES (qo = M D /M star ) AND GROOVE DEPTHS (A) CONSIDERED. 

Model q D A Qmin m Ju n (%i) fip {tip) 



1 


1 





1.21 


2 


1.21 (0.90) 


2.78 (2.69) 


2 


0.63 


0.90 


0.67 


2 


2.31 (2.14) 


3.41 (3.37) 


3 


0.32 


0.90 


1.31 


2 


1.43 (1.16) 


3.18 (3.05) 


4 


0.16 


0.90 


1.92 


2 


0.73 (0.65) 


3.04 (3.00) 


•5 


0.13 


0.90 


2.03 


2 


0.56 (0.50) 


3.02 (2.89) 


G 


0.08 


0.90 


2.24 


2 


(0.28 a - 6 ) 


(2.01) t 


7 


0.06 


0.90 


2.76 


2 


t 




12 


0.63 





1.30 


2 


0.84 (0.70) 


2.27 (2.22) 


13 


0.32 





1.52 


2 


0.37 (0.16, 0.31) 


2.03 (2.22) 


14 


0.16 





1.86 


2 


t 




15 


0.13 





2.08 


2 


t 




1G 


0.08 





2.26 


2 


t 




17 


0.06 





2.65 


2 


t 





TThis mode has two linear phases (a, b) and two "saturation" phases. The growth rates and pattern speed 
reported are measured from the hydro simulation. 

■'■These models do not show appreciable mode growth in either the linear or the fluid simulation. 
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Fig. 1. — Graphical depiction of the generalized eigenvalue problem. The square matrices represent the 
real and complex parts of the governing matrix. Solving the eigenvalue problem predicts the pattern speed, 
growth rate and radial dependence of the unstable mode. 
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Fig. 2. — Toomre Q — nc s /irGa profiles for the disk models considered (dashed lines: models with a groove, 
solid: models without a groove). Like colors correspond to equal qr>- 
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Fig. 3. — Surface overdensity S evolution for Models 1 and 2. 
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. — Top panel: normalized amplitudes of the m = 2 mode for disk models 1-3, 12, 13. Bottom panel: 
lized amplitudes of the m — 2 mode for disk models 4-6, 14-16. The slope of the dashed lines shows 
=>dicted prnwth rate of a. mode as resolved hv nnr linear mode code. 
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Fig. 5. — Evolution of azimuthally averaged density profile and density at the gap for models 5 (top panel) 
and 6 (bottom panel). 
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